Thermodynamically activated vortex-dipole formation in a two-dimensional 

Bose-Einstein condensate 
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Three distinct types of behaviour have recently been identified in the two-dimensional trapped 
bosonic gas, namely; a phase coherent Bose-Einstein condensate (BEC), a Berezinskii-Kosterlitz— 
Thouless-type (BKT) superfluid and normal gas phases in order of increasing temperature. In the 
BKT phase the system favours the formation of vortex-antivortex pairs, since the free energy is 
lowered by this topological defect. We provide a simple estimate of the free energy of a dilute Bose 
gas with and without such vortex dipole excitations and show how this varies with particle number 
and temperature. In this way we can estimate the temperature for cross-over from the coherent 
BEC to the (only) locally ordered BKT-like phase by identifying when vortex dipole excitations 
proliferate. Our results are in qualitative agreement with recent, numerically intensive, classical 
field simulations. 
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Since the experimental realisation of Bose-Einstein 
condensation (BEC) Qj, tremendous advances have been 
made in exploring the physics of ultracold atomic and 
molecular gases. It is now possible to simultaneously 
prepare multiple condensates and collide them 0, or 
even slice a condensate into numerous pieces using optical 
standing-waves |4|. Such an optical potential, in addition 
to the normal harmonic confinement, allows one to cre- 
ate a highly anisotropic trap and in this way it is possible 
to enter the quasi two- 0, IE IE or one-dimensional @ 
regimes. 

In lower dimensions the properties of quantum systems 
are cardinally altered. Unlike in the three-dimensional 
(3D) case, the different density of states means condensa- 
tion cannot occur either in two- (2D) or one-dimensional 
(ID) homogeneous ideal Bose gases in the thermody- 
namic limit at finite temperature [E@- 

The application, however, of a trapping potential in- 
troduces a low-momentum cut off which suppresses long- 
wavelength fluctuations and allows the possibility of BEC 
|lfj ITlT | below a certain critical temperature, which for 
the ideal 2D case we denote Tq D . 

In the homogeneous 2D system, although BEC is ab- 
sent at finite temperature, superfluidity can exist with 
algebraically decaying order fl^, El in an interacting 
but spatially uniform Bose gas. This is the Berezinskii- 
Kosterlitz- Thouless (BKT) phase in which the gas has 
local order but no global phase coherence, with no uni- 
form order parameter. A local order parameter can be 
defined and is known as the quasi-condensate 15] . 

A similar scenario persists in the presence of external 
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confinement. Three distinct regions can be identified in 
a two-dimensional trapped interacting bosonic system. 
At temperatures much below Tq D a phase coherent BEC 
appears to exist [Ifil, Il7j . At higher temperatures the 
coherence length decreases to length scales smaller than 
the system size [Tsj and bound vortex-antivortex pairs 
(vortex-dipoles) begin to form 19] indicating a BKT- 
like phase. Eventually, around Tq D , the vortex-dipoles 
unbind and the correlation length decays exponentially. 
This is the analogue of the BKT phase transition to the 
thermal gas |20ll21|. 

A phase diagram and approximate critical tempera- 
ture for the BEC to BKT cross-over were derived numer- 
ically in 0] assuming that the vortex-dipole is always 
at the centre of the trap, i.e. placed at the highest den- 
sity. Here, we use simple analytic arguments, without 
this restriction allowing the vortex-dipole to be placed 
anywhere inside the Thomas- Fermi radius, to obtain a 
very simple picture of the vortex dipole proliferation in 
the gas. Our approach obviates the need for any detailed 
numerical calculations and qualitatively agrees with the 
more sophisticated approaches, including the full classi- 
cal field simulations |20| . 

In our model, we minimize the Hclmholtz free-energy, 
F = U — TS, where U denotes the total internal en- 
ergy and S the entropy of the state, and the system has 
constant volume V and an absolute temperature T. We 
partition the gas into two subsystems, namely, a con- 
densate and a thermal part, each with its own internal 
energy and entropy contribution, (E c , E t ) and (S c and 
St), respectively. In the condensate part we allow for a 
vortex-dipole excitation, i.e. a bound vortex and antivor- 
tex pair. It suffices to consider only the effect of one pair 
in the condensate with the additional configurational en- 
tropy of the vortex dipole denoted as AS V . 
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The important contributions come from the conden- 
sate, with the thermal component only producing global 
energy shifts. Accordingly, E t and St are only roughly 
approximated. The free energy of the system can now be 
written as F = U - TS = (E c + E t ) - T(S C + AS V + S t ). 

The model system consists of N interacting bosonic 
particles of mass m confined by an external non-rotating 
harmonic potential characterized by the {lu±,uj z } trap 
frequencies. The actual values of the physical parame- 
ters we take correspond to the experiments of 0, 0] ■ In 
both experiments the radial and axial trapping frequen- 
cies differ by several orders of magnitude, suggesting that 
we use a quasi-2D formalism. This assumption is con- 
firmed by comparing the characteristic oscillator length 
and the s-wave scattering length Consequently, we 
assume that the wavefunction for the condensate can be 

I 

E c -N — + ^2y 



where V 2D = (d x ,d y ), and the potential term is 
cylindrically symmetric, (x s /a osc ,j_)(a; 2 +y 2 ). Interaction 
between particles is taken into account via the 3D s-wave 
scattering length, a s . 

The internal energy of the thermal part was approxi- 
mated by the energy of the trapped but non-interacting 
Bose-gas j2^| 

E -'»-*M !D (f)(jf)' <2) 

The corresponding entropy contribution St can be de- 
rived from Eq. J5J using the standard formulae of statis- 
tical physics. The validity of this approximation is also 
justified by experiment where the thermal cloud was 
not fully within the 2D regime. The total number of par- 
ticles was close to the theoretical threshold J2jj, above 
which the atoms squeeze out to the third dimension even 
if the condensate dimension in that direction is smaller 
than the healing length. 

Naturally, a more sophisticated treatment is possi- 
ble and in particular, the presence of the thermal cloud 
could significantly effect the energy of the vortex core 
|24| . To go beyond our simple treatment would require a 
numerically intensive, self-consistent, finite-temperature 
Hartree-Fock-Bogoliubov treatment [2^, l26l |27| which is 
not in the spirit of this exercise. Here we aim at the 
simplest possible estimate of the critical temperature for 
vortex dipole proliferation and we do not expect inaccu- 
racies in the vortex core energy at high temperatures to 
qualitatively effect our conclusions. This expectation is 
confirmed by comparison with the (numerical) classical 
field simulations of Simula and Blakie [20l |28| . 

In order to evaluate the free energy functional F, we 



factorized into two parts: $ = ip±(x,y)ij) z (z), where ip± 
describes the state in the xy-plane, while ip z is simply 
the ground state wavefunction of a harmonic oscillator 
in a trap characterized by u z . The wavefunction ip± we 
normalize to the condensate particle number Nq and ip z 
to unity. 

Moreover, we introduce an arbitrary spatial unit x s , 
rescaling the spatial variables and the wavefunction 
to be dimensionless. In the numerical calculation the 
spatial unit x s was taken to be equal to the oscillator 
length a osc ,j_ corresponding to the frequency luj_. 

Substituting this composite wavefunction into the 
usual Gross-Pitaevskii energy functional and evaluating 
all integrals with respect to the z direction, we obtain the 
following expression for the energy E c of the condensate 



dr, (1) 

I 

must also approximate the entropy, which is proportional 
to the number of physically available microstates. Only 
the ground state is available for particles in the conden- 
sate, therefore S c = 0. On the other hand, if a vortex- 
dipole is present in the condensate then its entropy can be 
estimated via a simple geometrical approach (see Fig. ^) 
such that AS V = &b ln(2x v /^(a; v )), where x v denotes 
the distance of the vortex-dipole from the centre. The 
healing length characterizes the length over which the 
density depression vanishes and the wave function ap- 
proaches its bulk value; it is defined by £ _1 / 2 (iz; v ) = 
4V2ir(a s /a OSCtZ )\(f> TF (x v )\ . 

It is now clear that only, E c and AS V , have spatial 
dependence and the remaining energy and entropy ex- 
pressions simply shift the frcc-cncrgy by a constant value 
at a given temperature. 
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FIG. 1: The dark ring represents the equivalent microstates 
for the vortex pair displaced at distance x v from the centre. 
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Finally, to evaluate E c in mean-field theory, one needs 
a wavefunction for the condensate containing a vortex- 
dipole excitation. The square modulus of the wave- 
function ipx without any excitation is chosen to have 
the Thomas-Fermi density profile on which we super- 
pose the vortex-dipole excitation in the form ijj±{v) = 
^>Tp(r)y-i(r)y>x( r )) where the vortex or ant i vortex wave- 
function is 



ip q (r) = tanh 



|r - r. 
£(r 9 ) 



91 \ e iq0(r,r q ) 



q = -!,+!. 



The two values of q correspond to a vortex and anti- 
vortex and r q — (x qi y q ) denotes the positions of the lo- 
calized excitations. The hyperbolic tangent ensures the 
proper asymptotes both in the vortex core and far from 
the vortex. The phase 9(r,r 9 ) measures the angle of r 
with respect to the x-axis. It gives a 2n phase differ- 
ence if one integrates along a path encircling the vortex 
core. We note here that the thermal cloud could infil- 
trate into the cores of vortex-di pole s and may enhance 
the expulsion of the condensate [24J which, as discussed 
previously, our model does not take fully into account. 
The vortex-dipoles appear at the edge of the condensate 
however, where the disturbance due to the thermal cloud 
is expected to be small. We therefore expect our simple 
model to be qualitatively accurate. 

Three different vortex-dipole configurations were ana- 
lyzed, namely a) = (x\,-yi), b) (x_i,0) = 
(-xi,0) and c) (a;_i,0) = (x x - £(xi),0). 

In the calculation, the total particle number N and the 
scaled temperature T/Tq D were varied. For every pair of 
N and T we determined the position x v — x\ of the vor- 
tex dipole to be the point where the free-energy reaches 
its minimum. These minima belonging to configuration 
a), are represented in Fig.^a). 

As expected, the position decreases as the tempera- 
ture approaches the Bose-Einstein condensation temper- 
ature Tq D , meaning that the minima of free energy are 
deeper inside the condensate. Simultaneously, however, 
the increasing particle number opposes the vortex-dipole 
formation well inside the condensate. Since vortex nu- 
cleation is not just a phase-pattern but also involves a 
density dip, more energy is required to create and main- 
tain a density disturbance in a larger condensate. 

In the T — > limit, one would anticipate that x v would 
approach i?TF- The small deviation from this expecta- 
tion seen in Fig.^a) is attributable to the numerical pro- 
cedure, since we require the vortex-dipole remain within 
the Thomas-Fermi radius. 

Theoretically speaking, statistical physics allows for 
the existence of a vortex-dipole excitation inside the con- 
densate phase with non-zero probability at every finite 
temperature, since the likelihood of nucleation is propor- 
tional to exp(— F/k-gT). However, there exists a well- 
defined temperature T c above which this probability be- 
comes comparable with unity. It is an important ques- 
tion how this temperature varies with respect to the total 
number of particles. Using the data depicted in Fig.J^a) 
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FIG. 2: (a) Scaled minima of free energy versus scaled 
temperature. Data are for configuration a). The horizon- 
tal dashed line is a reference used in the text, (b) The tem- 
perature required for a vortex-dipole to penetrate into the 
condensate to a distance of a function of total particle 

number N. The dashed line is a guide to the eye. 



one can infer the behaviour of T c (N) . Let us fix x v at an 
arbitrary value x^ (in our case 0.86-Rtf)- Calculating 
the intersections of this line with the curves belonging to 
various particle numbers [see the horizontal dashed line 
in Fig. Ufa)], one obtains (N^Ti) pairs. This T% repre- 
sents the 'critical' temperature for which a system with 
Ni particles has a minimizing vortex-dipole position at 
the given distance x^ from the centre. If we increase 
the temperature the vortex-pair can enter further. 

The results are represented in Fig. EJb) and show that 
the temperature finally approaches the condensate crit- 
ical temperature as the total number of particles grows. 
The slope of this graph naturally depends upon the value 
of x^ 1 but the tendency itself does not and these findings 
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are comparable with |Tgj . 

One can also conclude that it is rather difficult to ob- 
serve the transition from BEC to BKT state experimen- 
tally. At low ./V the condensate would not be dependably 
detectable after ballistic expansion, while at high N the 
vortex-dipolc formation in the bulk condensate phase re- 
quires a high temperature, which also results in a small 
condensate and large fluctuations. 

In conclusion, our simple model captures the physics 
of previous more complex analyses. There are two con- 
tributions to the free energy which depend on the posi- 
tion of the vortcx-dipole; the energy of the condensate 
and the entropy change due to the excitation. One an- 
ticipates that the energy should monotonically decline 
to the value of Thomas-Fermi energy as x v approaches 
i?TF, since the vortices deform the Thomas- Fermi density 
to an ever-decreasing extent. Similarly, the magnitude of 
the entropy contribution also decreases. This is because 
since the characteristic size of the vortices £ oc n -1 / 2 , 
the healing length rapidly increases. Although the area 
of the disk W(x v ) of available points in which to place 
the vortices grows linearly with radius, this does not re- 
sult in more states available for occupation due to the 
increasing size of the vortex-dipole. Consequently, these 
two competing geometrical effects are bound to lead to a 
minimum in the free energy F. 

A similar calculation has been performed in 29] and it 
was found that for a given angular momentum carried by 
the vortex dipole, the energy curve has a minimum inside 



the condensate. However, a direct numerical comparison 
with that result is not possible due to the differing ap- 
proximations used. 

Summarizing our result we can say that (i) the con- 
figurations of vortex dipoles are not energetically equiv- 
alent; configuration a) seems to be able to enter the bulk 
condensate more easily than the others as the tempera- 
ture increases, (ii) All types of configuration suggest that 
vortex-dipoles prefer nucleation at the edge of the con- 
densate, and they can penetrate further in the vicinity of 
the normal gas phase-transition (ss 0.8 — 0.9 Tq D ). This 
result is consistent with and supports the possibility of 
vortex-dipole nucleation in the experiment Q at ENS, 
Paris and is also consistent with their recent observation 
of the cross-over from thermal gas to BKT- like phase ■ 
(iii) With increasing particle number, the temperature T c 
characterising the vortex-dipole nucleation also becomes 
higher and asymptotically reaches the critical tempera- 
ture Tq D of BEC in the ideal gas. Qualitatively, the same 
tendency was found in pfl ] using the numerically inten- 
sive classical field simulation, in which dynamical effects 
were taken into account inherently. 

The nature of the change from BEC to BKT super- 
fluid, i.e. whether this is a slow continuous change or 
corresponds to a sharp cross-over, remains an open ques- 
tion and is the subject of further investigations. 

We express our thanks for financial support from the 
Marsden Fund and from the C. N. R. S. of France. 
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